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ABSTRACT 



Aims. To derive carbon isotopic ratios from optically thin tracers in the central regions of the starburst galaxies M 82 
and NGC 253. 

Methods. We present high sensitivity observations of CCH and two of its ^^C isotopologues, C^"^CH and ^"^CCH, as 
well as the optically thin emission from C^*0 and ^^C^*0. We assume the column density ratio between isotopologues 
is representative of the ^^C/^'^C isotopic ratio. 

Results. From CCH, lower limits to the ^^C/"C isotopic ratio of 138 in M 82, and 81 in NGC 253, are derived. Lower 
limits to the ^^C/^^C ratios from CO isotopologues support these. ^^C^^O is tentatively detected in NGC 253, which is 
the first reported detection in the extragalactic ISM. Based on these limits, we infer ratios of ^®0/^*0> 350 and > 300 
in M 82 and NGC 253, respectively, and ^^S/^*S> 16 in NGC 253. The derived CCH fractional abundances toward 
these galaxies of ^ 1.1 x 10~* are in good agreement with those of molecular clouds in the Galactic disk. 
Conclusions. Our lower limits to the ^^C/^^C ratio from CCH are a factor of 2 — 3 larger than previous limits. The 
results are discussed in the context of molecular and nucleo-chemical evolution. The large ^^C/^^C isotopic ratio of the 
molecular ISM in these starburst galaxies suggest that the gas has been recently accreted toward their nuclear regions. 
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1. Introduction 

Isotopic ratios in the interstellar medium (ISM) of star- 
burst (SB) galaxies provide important clues on their nucleo- 
chemical evolution. In particular, the ^^C/^'^C isotopic ratio 
is believed to be a good tracer of the chemical evolution of 
the Galaxy (see e.g. Audouze 1985) This ratio is understood 
to be a direct measurement of the primary to secondary nu- 
clear processing (Wilson & Rood 1994). While ^^C is a pri- 
mary product of stellar nucleosynthesis, quickly produced 
via He burning in massive stars that can be formed in 
first generation metal-po or stars. i s a secondary nuclear 
product from ^'^C seeds ( Mevedfiggl I Wilson k Matteuccil 
[l99l iWilson fc Ro od 199^ . 

So far, only lower limits to the ^^C/^^C abundance 
ratios could be determined for some starburst galaxies. 
A value of ^^C/^'^C ^ 40 was derived toward the nu- 
clei of the nearby galaxy NGC 253 based on observa- 
tions of CN, CS, and HNC (iHenkel et al.l[l99l and fur- 
ther supported by CO, HON, and HCO+ data on M 82 
and N GC 4945 (|Henkel fc Mauers bergei<1993HHenkel et al.l 
119941) . Additional CN data on M 82 and IC 342 resulted in 
lower limits to the ratio of > 40 and > 30 respectively 
(jHenkel et al.lll998D . The CO a nd "CO multi- transition 
non-LTE study toward M 82 by ()Mao et al.ll2000D excludes 
I2(~;^i3(n ].^(^jQg below 25 and point toward a ratio > 50. 
Apart from these multi-molecule studies toward the nearest 
SB galaxies, and limited by sensitivity, only the ^^CO/^^CO 
ratio is available for a small sample of nearby sources 



OYoung & Sanders' ^ToM; 'Sag e k Isbell [T99lt lAalto et al.l 
1991; Casoh ct al. 1992; Paglio ne et 311120011) . However, the 
lines from the species used so far can be se verely affected b y 
optically depth effects as estimated by IHenkel et al.l (|1993l ). 
Thus optically thin species are crucial to derive accurate 
isotopic ratios. 

Ethynyl (CCH) has the brightest line in the 2 mm 
spectral scans carr ied out toward NGC 253 and M 82 
(|Martm et al.l l2006l Aladro et al. in prep.). Its hyperfine 
structure, unresolved in these objects, can help to con- 
strain the opacity through the observed lineshape. Within 
the Galaxy, CCH emission appears tightly bound to mas- 
sive s tar formation at all evolutionary stages (jBeuther et al.l 

MM)- 

In this paper, we present observations of CCH and two 
of its ^'^C isotopologues, and provide a more stringent lower 
limit to the ^^C/^^C isotopic ratio in the nuclear starburst 
ISM of M 82 and NGC 253. Additionally we present CO 
observations in the ^^0 substitution and the double iso- 
topologue i^ciSQ. The derived ^^C/^^C ratio from CO, 
though not as stringent as that from CCH, supports the 
idea of larger isotopic ratios in these starbursts than previ- 
ously determined. 



2. Observations and Results 

We have observed the J = 2 — 1 group of hyperfine (hf) 
transitions of CCH (174.635 GHz), C^^CH (170.491 GHz), 
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and "CCH (168.274 GHz), and the J = 1 - transition 
of C^^^O (109.782 GHz) and "ci^O (104.711 GHz) toward 
M 82 and NGC 253. Observations were carried out with 
the IRAM 30m telescope (Pico Veleta, Spain). In the case 
of M 82, observations were aimed toward the north-eastern 
molecular lobe at the offset position (+13", +7.5") from 
the center (aj2ooo = 09h55™51!9, (5j2ooo = 69°40'47"1). 
NGC 253 observations were aimed at the central posi- 
tion (aj2ooo = 00h47™33'!3, (5,72000 = -25°17'23"15) for 
C^®0 while the ^•^C^^O observations were slightly offset 
at (+3", —4"). The beam widths at these frequencies were 
- 14" and - 23", for CCH and CO lines, respectively. We 
used the wobbler switched observing mode with a symmet- 
rical beam throw of 240" in azimuth and a wobbling fre- 
quency of 0.5 Hz. As spectrometers we used the 256 x 4 MHz 
fiher banks and WILMA autocorrelator (2 MHz) for CCH 
and CO, respectively. The CCH data for NGC 253 were ex - 
tracted from the 2 mm line survey bv lMartm et al.l (|2006D . 
Figs. [1] shows the observed spectra. 

The hf structure of CCH is unresolved. We fitted a 
comb of Gaussian profiles at the positions of the hf lines 
and intensity ratios fix ed by their spectroscopic parame- 
ters (|Muller et al.ll200ll ). shown as vertical lines in Fig. [T] 
In M 82, the wing of emission at low velocities in CCH is 
a contribution of both the fainter hf components and the 
emission from the nuclear region, as also seen in the C^*0 
profile. Though two velocity components could be fit to 
the CCH profile in NGC 253 (jMartm et al.ll2006l) . for the 
purposes of this paper we just fit a single velocity com- 
ponent. We assumed optically thin emission to fit the line 
profiles. We show the resulting fits as grey lines in Fig [1] 
while the fitting parameters are shown in Table [T] Neither 
C^^CH nor ^^CCH emission were detected toward M 82 
and NGC 253. No line contamination by other species is 
expected at the frequencies of the isotopologues. We have 
estimated upper limits to the integrated intensities of the 
overall line profile of these isotopolog ues, calculated from 
the corresponding spectral parameters ()Miiller et al.ll200"l[) . 
We assumed a peak intensity of the overall profile at a 3<t 
level (for 187 and 224 kms~^ wide channels, corresponding 
to that from single Gaussian fit to the CCH features) for 
M 82 and NGC 253, respectively. 

Line emission from ^'^C^^O is not detected toward 
M 82 and only tentatively toward NGC 253. This tenta- 
tive ^■^C^^O detection, the first in the extragalactic ISM, 
appears blended to the H2CS 3i,2 — '^1,1 transition at 
104.616 GHz. The derived intensit y of ^ 1 mK suppo rts 
the tentative detection reported bv iMartin "eTall (pOOSh . 

With the fitted integrated intensities we have estimated 
the beam averaged column densities for each species un- 
der the local thermodynamic equilibrium assumption. An 
excitation temperature of Tqx =10 K was assumed for the 
calculations. The assumed temperature , being the same for 
all isotopologues, has no effect on the derived column den- 
sity ratios. As reported from CO observations, the differ- 
ent isotopologues might be tracing different gas component 
(jWal l & Jaffc 1990; Aaho et al. 1994; Meier et al. 2000), 
and therefore the similar excitation temperature assump- 
tion would not be fulfilled. As a test case, if the ^'^C isotopo- 
logue was emitted from a gas component with T^x '^50 K, 
our assumption would lead to the ^^C/^'^C ratio being over- 
estimated by a factor of ^ 1.7. Multi transition observations 
of CCH and its ^'^C isotopologues would be needed to accu- 
rately establish the magnitude of such effect. On the other 
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Fig. 1. Observed spectra of the J = 2 — 1 transitions of 
CCH, C^^CH, and ^^CCH, and the J = 1 - transitions of 
C^^O and ^^C^^O toward M 82 and NGC 253. The intensity 
of the C^^O spectrum in M 82 has been multiplied by a fac- 
tor of 3 for comparison and visualization purposes. The fit 
to the CCH and C^^O profiles are shown with grey lines. For 
the undetected isotopologues we overlay the expected spec- 
tra assuming a ^^C/^^C isotopic ratio of 40 (dotted lines), 
as well as the upper limit profile derived in this work (grey 
lines). Vertical lines represent the relative intensities of the 
hyperfine components used to fit the detected CCH profiles, 
and to calculate the upper limit to the intensities for the 
isotopologues. ^^C^^O is tentatively detected blended to a 
H2CS feature toward NGC 253. Spectra are presented at 
the original 4 and 2 MHz (~ 7 and 5km s~^) resolution for 
CCH and C^^O, respectively. C^^CH and ^^CCH data were 
smoothed to 16 MHz (- 28kms-i) and ^^C^^O to 8 MHz 
(- 23kms-i). 



hand, we do not expect it to be significant in the optically 
thin CO isotopologues. The derived CCH and CO beam 
averaged column densities and upper limits are shown in 
Table m 

The derived lower limits to the ^^C/^^C ratios are of 
> 110 and > 56 for M 82 and NGC 253, respectively. We 
can further constraint these limits by averaging the spectra 
of both C^^CH and ^^CCH. The resulting limit to the col- 
umn densities raises the limits to the isotopic ratios up to 
i2Qy'i3Q ^ ]^3g g^j^^j ^ g]^^ respectively. Using CO we obtain 

ratios of > 56 and ^ 60, respectively. 
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Table 1. Observed lines parameters, and derived column densities and ^^C/^'^C ratios 
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Notes. This value refers to the brightest component in the group. la rms at the resolution shown in Fig. [T] ''^' Beam 
averaged column density assuming a Tex = 10 K. '''' 3a limit to the integrated intensity '"^^ Tentative detection blended to an H2CS 
transition. Velocity and width parameters were fixed to fit both lines simultaneously. 



3. ^^C/ "^C isotopic ratio from CCH and CO 

Estimating the interstellar isotopic ratio from the derived 
molecular abundances one has to discuss two main draw- 
backs, opacity and fractionation effects. 

Lines of ^^C isotopologues of CO and abundant species 
such as HCN or HCO+, are likely to be optically thick. If 
our derived lower limits to the carbon isotopic ratio ap- 
plies to CO, it would imply opacities of ti2co 5 and 
^ 9 for NGC 253 and M 82, respectively, as derived from 
the observed CO to ^'^CO integrated intensity ratios by 
ISage fc Isbelll (|1991[ 1. Higher opacities will lead to a de- 
crease of the measured ^'^C/^^C ratio from CO and possibly 
other abundant species. Therefore, molecules with opaque 
lines poorly constrain the ^^C/^^G ratios. From the hf fit 
to the line profile we can exclude large opacity effects af- 
fecting the CCH lines. Moreover, small opacity effects on 
the CCH observed lines would result in an increase of our 
lower limits to the ^^G/^^G isotopic ratios. 

Fractionation (IWatson et all 119761: iLanger et al] 11984 
IWilson fc Roodlll994 ) might produce an enhancement of 
-"^•^CO in the outer layers of molecular clouds. Thus the ob- 
served CO would be tracing regions with lower ^^CO/^^^CO 
ratio than the actual ^ ^C/^'^C isotopic ra tio. The CCH ob- 
servations reported by ISakai et aD (|201Cl[ l toward the dark 
cloud TMC-1 and the star forming region L1527, also show 
a ^^C/^'^C ratio a factor of 2 — 4 larger than the inter- 
stellar ratio. However, the chemical fractionation claimed 
in their work is mostly effective at very low temperatures 
(Tkin ^ 10 K) while isotopic exchange re actions rates tend 
to ba lance out at higher temperatures ([Woods fc Willac"vl 
l2009f ). This fractionation is also expected to affect species 
other tha n CO, such as HCO"*" an d HCN. Moreover, the 
models bv lWoods fc Willacvl (|2009[ ) show how species such 
as HCN and HNC should result in even larger ratios than 
those from CCH, which is not observed in galaxies. Thus, 
chemical fractionation cannot explain the large CCH ratio 
observed in the warm ISM in the central regions of galax- 
ies. ^^C/^^C isotopic ratios of 20 — 25 derived from CCH in 
the GC (Armijos et al. in prep.), similar to those derived 
from other molecular species, support the hypothesis that 
fractionation does not play an important role. 



We conclude that ours limits for the ^^C/^'^C isotopic 
ratios derived from CCH isotopologues are representative 
of the ISM in the observed SB galaxies. 



4. A large ^^C/^-^C isotopic ratio in starbursts 

Our results from CCH show that the ^^C/^'^C isotopic ra- 
tio in the SB environment is much higher than the value 
of ~ 20 measured in the Galactic Center (GC) region 
(jWilson fc RoodI |1994[ ). Moreover, our limits are clearly 
higher than p revious limits for SB galaxies of ^ 40 — 50 
(|Henkel et al . 1993; H enkel fc Mauersbergedll993D . We can 
exclude a ratio of 40 or lower, since in that case we would 
have detected the ^^CCH emission in both galaxies and 
that of C"CH in M 82. If the ^^G/^^C ratios derived from 
CCH are representative of the overall '^'^G/'^^G, it would 
imply this ratio to be at least a factor of 2 — 3 larger than 
previously reported in the starburst ISM. Though lower, 
the limits to the ratios we derive from CO isotopologues 
also suggest larger ^^C/^^C isotopic ratios than previously 
measured. 

Single dish ^^CO/^^CO ratios have been measured as 
high as 55 and 44 toward N GC 4195 and NGC 6240, re- 
spectively (jCasoli et al.lll992l ). However the latter was dis- 
proved b y further mapping resulting in a ^^CO/^'^CO ratio 
of ~ 20 |Aalto fc Hiittemeisteill200tf ). In Arp 299 (consist- 
ing of a merger of IC 694 and NGC 3690), a ratio of 60 
was found toward the nucleus of IC 694, and ~ 10 in the 
surrounding disk (jAalto et al.lll999t ). with the ratio rang- 
ing over ^ 20 — 35 at single dish resolution. Recently, a 
i^CO/i^co J = 3 - 2 ratio of > 40 was reported toward 
the Clover leaf quasar at z ~ 2.6 (jHenkel et al.ll2010l ) 



5. Revisited carbon isotopic ratio implications 

5.1. The oxygen i^O/^O ratio 

The uncertainty in the measured ^'^G/^^G ratio directly af- 
fects the determination of the ^^O/^^O ratio based o n mea- 
surements of "CO and C^O. For NGC 253, Harriso n et all 
()1999D presented ^'^CO and C^^O measurements. If we use 
our revised lower ^^G/^^C limit of > 81, their "CO/C^^O 
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intensity ratio would correspond to an abundance ratio of 
16q/1 8q^ 300 Toward M 82, using the ^^CO observations 
from (|Mao et al.ll2000f ). observed at the same position, we 
derive a ^^0/^®0> 350. Both Umits are now closer to the 
ratio of ~ 3 30 found within the 4 k pc molecular ring of the 
Milky Way (iWilson fc RoodllTool . 



5.2. The sulfur ^^Sf-^S ratio 

The ^^S/^'*S ratio based on measurements ofJfCS_aiid CffS 
will also be affected, as that derived by Ma rtin et al.l (|2005[ ) 
for NGC 253. Their value of the ^^S/^'^S ratio of 8 ± 2 
would be doubled up to a value > 16, closer to the valu e 
of ^ 24 measured in the Galactic disk (IChin et al. 1996fl. 
The ratio of 32s/34s- 13.5 toward NGC 4945 ("Wan g et all 
[2003) might also be underestimated as it was based on the 
CS/C'^^'S ratio assuming no significant saturation of the CS 
lines. 



5.3. H2 column density determination 

Our large lower limits to the ^^G/^^C ratio impacts the H2 
column density {NH2) derived from the optically thin CO 
isotopologues. These isotopologues, less affected by optical 
depth, are expected to be better suited to trace the overall 
molecular material than the main isotopologue. However, 
the determination of the will depend on the isotopic ra- 
tios. If we apply our new limits to the oxygen isotopic ratio 
to estimate the from C ^^O, and assum ing a conversion 
factor CO/H2 ~ 8.5 x IQ-^ (jFrerking et aWmS'i) . we derive 
beam averaged hmits of 7Vh2 > 2.6 x 10^^ cm"^ for M 82 and 
7Vh2 > 1.2 X 10^3 cm-2 for NGC 253. We then infer CCH 
fractional abundances of < 2 x 10~* and < 0.6 x 10~*, 
respectively, in good agreement with the Galactic abun- 
danc es found in prestellar cores (1 x 10~^, Padovani et al.l 



I2009D . dark clouds (< 6 x IQ-MW ootten e t al.lll980D . dif- 
fuse clouds (3 x 10-Ml ucas fc Lisz t 2000) and about two 
orders of mag nitude above the abun dances in hot cores 
(2 - 8 X 10-i°. lNummelin et al.l[20(B . 

Moreover, assuming our derived ratio of ^^O/^^O^ 300 
a better agreement is found between the gas mass estimate 
from 1.3 mm dust observations (Krueel et al, .1990) and 
that from C^^O (see Table 2 in iMauersberger et al.lll996l 
where a ratio of ^^0/^*0~ 150 was used). This agreement 
further supports our result. 

We can estimate a conversion factor of Xco — 



conversion 



7V(H2)/CO - 1.3xl0^"cm-^(Kkms-^ ),using the CO 1-0 
integ rated intensity of 920Kkms~^ (jMauersberger et al.l 
Il996|) . This conversion factor, significantly lower by a factor 
of ~ 2 than the "standar d" value within the Galaxy of 2 — 
^fKkms"^) (|Solomon et al.lll987HStrong et all 
factor of 3 larger than the values proposed 



3x10^° cm 



I1988D . IS a 

for the starburst ISM of 0.3 



0.4 X 102°cm -2(Kkms"^) 



(jMauersberger et al.l[T996HMartfn et al.ll2010[ ) 



6. Chemical evolution of starburst galaxies 

After several cycles of star formation, nucleosynthesis will 
enrich the interstellar medium with processed material 
which could be "dated" by using the ^^C/^'^C isotopic ra- 
tio. Thus, the ^^C/^^C ratio is exp ected to decrease with 
time. The solar system ratio of 89 (jWilson fc Roodlll99l ) 
is thought to be representative of the local ISM when the 



Solar System formed, 4.6 billion years ago. Although the 
i^CO/^^CO ratio toward both M 82 and NGC 253 is ob- 
served to be a factor of ~ 1.7 larger within the nuclear 
region than in the outer disk, this ratio is unlikely to be 
representative of the ^"^C/^^C ratio due to opacity effects 
and/or va riation of the CO/H2 conversion factor, as re- 
ported bv iPaglione et al.l (|200l[ ). On the contrary, a pos- 
itive gradient towards the central region is expected, due 
to the molecular material being further processed by stars, 
and similar to wh at is observed with optically thin tracers 
in th e Milky Way (|Langer fc Penziadll990H Wilson fc RoodI 
'1994). It is remarkable that, taking into account the ra- 
tios derived in this work, the ^^C/^^C ratio in SB galaxies 
is lar ger than in the Galac tic center, the solar neighbor- 
hood i|wilionXRoo3[l99l and even in the low metallicity 
ISM of the Large Magellanic Cloud (LMC; ^^C/^^C= 49, 
IWang et al.ll2009D . Similarly, the oxygen and sulfur isotopic 
ratios measured do not evidence an enrichment of ^^0 and 
^''S in the starburst IS M with respect to the values mea- 
sured in the local ISM (jWilson fc RoodI I1994D . Moreover, 
our lower limit to ^^O/^^O is slightly ab ove the value of 
250 m easured toward the Galactic Center (jWilson fc RoodI 
I1994D . and the ^'^S/^^S would be at least a factor of 2 
above the expected ratio within t he central 3 kpc o f the 
Galaxy from the gradient found by IChin et al] (jl996|) . The 
16q^18q g^j^^ 32g^34g ^^.^ significantly lower and similar, re- 
spectively, to those measured in the LMC (^^0/^^0=2000, 
32s/34s^ 15. IWang et J1l2009h . 

An enhanced of star formation was likely triggered 
across the disk of M 82 as a result of a close encounter with 
M 81 about 220 Myr ago. However, the ongoing nuclear 
starbursts started ^ 15 Myr ago, with nuclear star clus- 
ters dated between 7 and 15 Myr (jKonstantopoulos et al.l 
,200!) . This dating is in agreement with the two short du- 
ration bursts of star formation i n the nuclear region of 
M 82 , with peaks at 5 and 10 Myr (jForster Schreiber et al.l 
I2003D . The nucleus of NGC 253 appear to host an 
e ven younger starburst with a median age of ~ 6 Myr 
(Ferndndez-Ontiver os et al.ll2009D . With a average lifetime 
of 10 Myr, massive star clusters formed in the starburst ap- 
pear not to have efficiently enriched the ISM in their nuclear 
region. Furthermore, the enriched material could have been 
banished from the nuclear region via the starburst driven 
super- winds in both galaxies (Hcckman ct al. 199(D. 

One possible interpretation would be that in very young 
SBs one should find large ^'^C/^^G ratios due to fast evo- 
lution of very massive stars, leading to a overproduc- 
tion of ^^C relative to ^^C, synthesized in intermediate 
mass stars (Henkel fc Mauersber ger 1993). However, such 
fast evolution might also result in an overproduction of 
^^O ([Henkel fc Mauersberger|[T993D which is not observed. 
Moreover, if the time scale of the main burst in the nu- 
cleus of M 82 is less than 15 Myr ((Fdrster Schreiber et al.l 
i2003 : Konstantopoulos et al.ll2009D . then most of the 
enriched material would be in the hot outflow with con- 
tribution to the molecular gas. The time scale is rather 
short to have the ^^C enriched material well mixed over 
the ~300 pc of our beam size, in both M 82 and NGC 253. 
It would mean that the enhanced ^'^C molecular gas is 
confined to shells around the stellar clusters with a huge 
12C/"C ratio. 

Thus, our measured isotopic ratios suggest that the bulk 
of the mass of molecular material within the nuclear regions 
of M 82 and NGC 253 consists of little or unprocessed gas 
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by massive stars. This suggests that the ISM which is sup- 
porting the starbursts in M 82 and NGC 253 has been re- 
cently funneled toward their centers from the outer regions 
or accreted by these galaxies. The evolution of these iso- 
topic ratios as a function of meta llicity is not clear fr om the 
available observations. However, IWang et al.l ()2009[ ) found 
the ^^O/^^O ratio to be a good tracer of metallicity. The 
lower limits to the ratio found toward M 82 and NGC 253, 
larger than the value toward the Galactic Center, is consis- 
tent with that found within the 4kpc molecular ring and 
even that within the local ISM. This would support the idea 
of inflowing material from the outer regions of these galax- 
ies toward their centers conforming the gas concentration 
unleashing their starburst events, where enriched material 
from previous star formation events would contribute little 
to the bulk of the starburst molecular material. 
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